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Abstract. A variety of terrestrial planets with different physical parameters and exotic atmospheres might plausibly 
exist outside our Solar System, waiting to be detected by the next generation of space-exploration missions. Some of 
these planets might transit their parent star. We present here the first study of atmospheric signatures of transiting 
Earth-size exoplanets. We focus on a limited number of significant examples, for which we discuss the detectability 
of some of the possible molecules present in their atmospheres, such as water (H2O), carbon dioxide (CO2), ozone 
(O3) or molecular oxygen (O2). To this purpose, we developed a model to simulate transmission spectra of 
Earth-size exoplanets from the ultraviolet (UV) to the near infrared (NIR). According to our calculations, the 
signatures of planetary atmospheres represent an absorption of a few parts-per-million (ppm) in the stellar flux. 
The atmospheres of a few Earth-like planets can be detected with a 30-40 m telescope. The detection of the 
extensive atmospheres of tens of small satellites of giant exoplanets and hundreds of hypothetical ocean-planets 
can be achieved with 20-30 m and 10-20 m instruments, respectively, provided all these planets are frequent and 
they are efficiently surveyed. We also found that planets around K stars are favored, mainly because these stars 
are more numerous and they are smaller compared to G or F stars. While not addressed in this study, limitations 
might come from the stellar photometric micro-variability. 
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^ I 1. Introduction as they were transiting their parent star, producing a pho- 

r* ■ tometric occultation. The last transiting planet identified 

t > ■ The Earth is the only known example of a life-hosting is a Saturn _ mass planet orbiting HD 149 026, a bright 

k> , world, even though terrestrial exoplanets have been y = g 15 GOlv star ( Sato et aL 20 05). The first discov- 

■ searched for since the discovery of the first Earth-mass ered transiting giant exoplanet, HD 209458b (Henry et al. 

. CO ; exoplanets by Wolszczan & Frail (1992). However, planets 2 000; Charbonneau et al. 2000; Mazeh et al. 2000), is the 

similar to the Earth, Venus or Mars, in size, density or or- object of intense investigations dedicated to characterizing 
bital parameters are still beyond the reach of the present nQ ^ a t m osphere 

capabilities for planet detection around normal stars. 

TT Probing planetary atmospheres by stellar occultations 

Until now, mostly giant exoplanets have been discov- . — ,. ,, , , r , , P , , . 

. _ ... . . , ... is an effective method used for a lot of planets and their 

ered. Remarkable progresses have been made recently with , .... . „ . „ „ T , , „. , 

. .. r , • r-,A, „, -n i satellites m the Solar System, from Venus to Charon (see, 

the discovery ol planets m the mass range ol 14 to 21 Earth , . . lr .„„ % „,.., , 

/■, , , „, t, r a , , r, e.g., Elliot & Olkm 1996). With this technique, we can ob- 

masses (14 to 21 M ffl , see McArthur et al. 2004; Santos et ' , , • . A 

. _ _ _ ,\ . . serve the thm atmospheric ring surrounding the optically 

al. 2004), and most recently a ~7.5 M ffi planet orbiting . , ,. . , . , .. . T . . 

o^r. , ttt 1 , , thick disk of the planet: the limb. In the case of giant ex- 

G J 876 Rivera et al. 2005). We may speculate, then, that , , , . , . , /1 nM 

V ... . , , _ , ' . . . oplanets, though, the star is only partially occulted (1.6% 

smaller planets with sizes down to that ol the Earth might „ ... , tT mnnnirDn ™ , <• 

, , . . _ . . . , °. for the transiting planet HD 209 458b). The spectrum of 

be observed m a near luture. Among the 161 planets de- ,. , , . ' . , , , , C1 , , , , , 

the star light transmitted and filtered by the lower and 

tected so far, eight have been discovered or re-discovered „., . , , , „ 

truck giant exopfanct atmosphere consequently presents 

extremely weak absorption features (from 10~ 3 to 10~ 
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2001). 
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atwww.ciw.edu/boss/IAU/div3/wgesp. detected: Charbonneau et al. (2002) measured the 



at vo.obspm.fr/exoplanetes/encyclo/encycl.html. See also 
the web page of the IAU Working Group on Extrasolar Planets Despite the difficulties, such dim signatures were 
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lower atmosphere of HD 209 458b as they detected a 
(2.32 ± 0.57) • 10~ 4 photometric diminution in the sodium 
doublet line of the parent star at 589.3 nm. However its up- 
per atmosphere, which extends up to several planet radii, 
shows even larger signatures. Vidal-Madjar et al. (2003, 
2004) observed a 15 ± 4% absorption in the Lyman a 
(Ly Q ) emission line of HD 209 458 at 121.57 nm as well 
as absorptions from atomic carbon (7.5 ± 3.5%) and oxy- 
gen (13 ± 4.5%) in the upper atmosphere. In this work, 
we will discuss the possibility to detect and to character- 
ize the lower atmospheres of exopianets using signatures 
comparable in origin to the one detected by Charbonneau 
et al. (2002). 

The idea is to extend the use of transmission spec- 
troscopy to hypothetical Earth-size planets. We estimate 
that these exopianets present at least two orders of mag- 
nitude less signal than gaseous giants, as the transit of 
the planet itself would have a dimming of ^10~ 5 (the 
transit depth, AF/F, where F is the stellar flux, can 
be expressed as (Rp/R*) 2 , with Rp and R± standing 
for the radii of the planet and the star, respectively). 
The atmospheres of Earth-size exopianets should span 
over ~ 100- km height without considering potential upper 
atmospheres. Depending on their transparency - which 
would give an equivalent optically thick layer of ~10 km - 
the expected occultations caused by atmospheric absorp- 
tions should be ~10" 7 to ~10~ 6 . 

Earth-size planets are probably the most challenging 
objects to detect with transmission spectroscopy. The or- 
ders of magnitude given above, in fact, raise many ques- 
tions: is it realistic to seek for possible features that dim, 
with an instrumentation that might or might not be avail- 
able in a near future? What are the strongest signatures 
we should expect? What kind of planet could be the best 
candidate to look at? 

We have developed a one-dimensional model of trans- 
mission at the limb to give quantitative answers to these 
questions. Since we use the stellar light to explore the plan- 
etary atmospheres, we chose to focus on the wavelength 
range where the largest number of photons is available, 

1. e., between 200 and 2 000 nm. The model is described in 
Sect. The detectability of the selected atmospheres de- 
pends on the signal-to- noise ratio (S/N) achievable with 
a space telescope spectrograph. The constraints on ide- 
alized observations and the method we used to calculate 
their S/N, are described in Sect. OH Finally, the results for 
the specified cases are given and discussed in Sect. 0] 

2. Model description 

2.1. Geometric description of the model 

The general geometry of a transiting system is described 
by Brown (2001). In the present work we consider a non- 
transient occultation for the 'in transit' phase, with a null 
phase configuration (configuration 2 in Brown's Fig. 1), 
that is, the planet is centered in the line of sight with 
respect to the star. This configuration both maximizes the 




atmosphere 



Fig. 1. Sketch of the transmission of the stellar light 
through the planetary limb. The planet itself, i.e. the 
'solid' disk (in grey) is optically thick at all wavelengths. 
The quantity dl is the elemental length along the line of 
sight. In the calculation, we prefer to use the height h in- 
stead of I. The scale in the figure has been distorted for 
clarity. 

area of the atmosphere that is filtering the stellar light and 
minimizes any effects linked to the stellar limb darkening 
(Seager & Sasselov 2000). 

The stellar light is filtered through the atmospheric 
limb of the planet, as sketched in Fig.^ In the following we 
detail the integration of the atmospheric opacity along a 
stellar light path (or cord) through the limb of the planet. 

2.1.1. Opacity along the line of sight 

We calculate the total opacity of the model atmosphere, 
t\, along a cord, parallel to the line of sight, as the sum 
of the opacity of each species i present in the atmosphere, 
t\ = J^i t a,i- We can calculate the opacity along the cord 
as a function of its impact parameter, b: 

r+oc 

T X ,i(b) = 2 / A x , iPi (h)dl, (1) 
Jo 

where is the absorption coefficient for the species i at 
the wavelength A, expressed in cm 2 g _1 , and Pi(h) is the 
mass density in g cm" 3 of the species i at an altitude h in 
the atmosphere. 

Now, re-expressing Eq. ^ as a function of the height 
z = h + Rp, with Rp being the planet radius, we obtain: 

rx,i(b) = A KlPl (z - R P ) -0L= , (2) 

where £> max is the height of the higher atmospheric level we 
are considering. The method to estimate 6 max is presented 
in Sect. l2~2~31 

2.1.2. Spectrum ratio 

Consider the stellar flux received by the observer during 
the planetary transit to be F lni and the flux received when 
the planet is not occulting the star to be F out . Brown 
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(2001) defined 5i to be the ratio between those two quan- 
tities, and 3?' (the so-called spectrum ratio) as 3?' = 5R— 1. 
Here, 3?' is the sum of two distinct types of occultations: 

— The occultation by the 'solid' surface of the planet, 
optically thick at all wavelengths. Projected along the 
line of sight, this is a disk of radius Rp and the occul- 
tation depth is simply (i?p/i?*) 2 . 

— The wavelength-dependent occultation by the thin ring 
of gaseous components that surrounds the planetary 
disk, which can be expressed as T,\/(ttRI). The area, 

is the atmospheric equivalent surface of absorption 
and may be calculated as: 







/ 2ttM6 


l_e-^W' 


'Rp 





The resulting spectrum ratio is: 

#(A) = -^±^. (4) 
Note that $t? < 0. 

2.2. Description of the atmospheric profiles 

Along a single cord, stellar photons are crossing several 
levels of the spherically stratified atmosphere. We gener- 
ate an atmospheric model using the vertical profiles from 
Tinetti et al. (2005a, 2005b) and Fishbein et al. (2003) 
for the Earth and from the Venus International Reference 
Atmosphere (VIRA, Kliore et al. 1985) for Venus. These 
atmospheric data include the profiles of pressure, p, tem- 
perature, T, and various mixing ratios, Y. The atmo- 
spheres are initially sampled in 50 levels, ranging from the 
ground level to an altitude of about 80 km for the Earth 
and about 50 km for Venus. Both profiles stop below the 
homopause, so we assume hydrostatic equilibrium for the 
vertical pressure gradient. 

A useful quantity to describe atmospheres in hydro- 
static equilibrium is the scale height, H, i.e. the height 
above which the pressure decreases by a factor e. The 
scale height explicitly depends on the temperature, as 
H = UMaT / (p,g), where k and Na are the Boltzmann's 
and Avogadro's constants while p is the mean molar mass 
of the atmospheric gas. Since g is the acceleration due to 
gravity, H also implicitly depends on the radius and the 
density of the planet 2 . Consequently, less dense objects 
are likely to have more extensive atmospheres, hence they 
are easier to detect (Brown 2001). 

Density and size of planets are therefore key param- 
eters for the present work. In order to estimate their in- 
fluence, we test a set of different planetary types ranging 
from the Titan-like giant planet's satellite (pp « 2 gcm~ 3 , 
Rp ss 0.5 Earth radius - 0.5 R®) to the 'super-Earth' ob- 
ject (pp w 6 gem" 3 , Rp w 2 Rq). For the physical prop- 
erties of plausible, theoretically predicted planets such as 

2 To avoid confusion between the density of the atmosphere 
and the mean density of the planet, the latter is denoted p P 



a 'super-Earth', we use the mass-radius relation model 
from Dubois (2004) and from Sotin et al. (2005). Our at- 
mospheric model allows the re-scaling of vertical profiles 
depending on the acceleration due to gravity of the planet 
and the atmospheric pressure at the reference level. 

2.2.1. Molecular composition of the atmosphere 

Our simplified atmospheric profiles contain only the 
species that may produce interesting spectral signatures 
in the chosen wavelength range (0.2 to 2 pva), viz., wa- 
ter vapor (H2O), carbon dioxide (CO2), ozone (O3) and 
molecular oxygen (O2). Molecular nitrogen (N2) has also 
been considered, though lacking marked electronic transi- 
tions from the UV to the near IR. Nevertheless, it is a ma- 
jor species in Earth's atmosphere and it has a detectable 
signature via Rayleigh scattering at short wavelengths. 

We consider three types of atmospheres: (A) N2/O2- 
rich, (B) C0 2 -rich and (C) N 2 /H 2 0-rich cases. The first 
two types can be associated with existing planetary at- 
mospheres, respectively Earth and Venus. The last type 
(C) could correspond to the atmosphere of an Earth-mass 
volatile-rich planet such as an 'ocean-planet' described by 
Leger et al. (2004). The basis for building a 'toy model' of 
an H2 0-rich atmosphere are found in Leger et al. (2004) 
and Ehrenreich et al. (2005b, see Sect. 12.2. 5|) . 

Vertical gradients in the chemical composition and 
temperature of each of these atmospheres are plotted 
in Fig. [21 (N 2 /0 2 -rich), Fig. (C0 2 -rich) and Fig. H 
(N2/H2 0-rich). Table ^ summarizes the mean chemical 
compositions of these model atmospheres. 

2.2.2. Temperature profiles 

As mentioned above, we use Earth and Venus vertical tem- 
perature profiles as prototype for N 2 /02-rich and CO2- 
rich atmospheres (see Sect. I2.4|l . Moreover we assume an 
isothermal profile in the thermosphere, instead of the real 
one. This is an arbitrary, but conservative choice, since 
the temperature should on the contrary rise in the ther- 
mosphere enhancing the atmosphere's detectability (see 
Sect. 1323. 

2.2.3. Upper limit of the atmosphere 

We set the profiles to extend up to a critical height 6 ma x 
from the centre of the planet, or /i max from the surface 
(frmax — Rp + ft max ). This limit corresponds to the alti- 
tude above which the molecular species we considered 
(H 2 0, O3, CO2, O2) are likely to be destroyed or mod- 
ified either by photo-dissociating or ionizing radiations, 
such as Ly a or extreme-UV (EUV). 

Therefore, the critical height corresponds to the 
mesopause on Earth (at « 85 km). The column den- 
sity of the terrestrial atmosphere above that altitude, 
A/>85km 7 is sufficient to absorb all Ly Q flux. In fact, as 
the number density of the atmospheric gas, n(h), de- 
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Type fi (gmo-r 1 ) Fn 2 (%) Yh 2 o (%) Y C p 2 (%) Yp 2 (%) Yp 3 (%) Used for models 
N 2 /0 2 -rich 28^8 78 03 O03 21 < lO^* Al, A2, A3 

C0 2 -rich 43.3 4 3 • 1(T 4 95 Bl, B2, B3 

N2/H 2 0-rich 28.7 80 10 10 CI, C2, C3 

Table 1. Mean volume mixing ratio of atmospheric absorbers for the different types of model atmospheres considered. 
(*) Ozone is only present in model Al. 
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Fig. 2. Atmospheric profiles, Al. The plot shows the total 
number density profile (thin solid line) of the atmosphere 
in cm -3 , and that of the five species included in our model, 
namely, N2 (dotted line), O2 (dash-dot-dot-dotted line), 
H2O (dashed line), CO2 (long-dashed line) and O3 (dash- 
dotted line). Temperature (thick line up to 80 km) and 
mixing ratios of the different species are those of Earth. 
Temperature is assumed to be constant above that height. 
The thickest horizontal line shows the position of the cloud 
layer. 

creases exponentially with height, we can simply consider 
A/">85 km °C 7^85 km • #85 km, where 

^85 km 

and _ff 85 

km are 

the density and the scale height of the terrestrial atmo- 
sphere at 85 km, respectively. 

Similarly, we set the upper limit of a given atmo- 
sphere, h max , to the altitude below which the photo- 
dissociating photons are absorbed. We assume that 
hm&x is the altitude where the column density equals 
that of the terrestrial atmosphere at 85 km, that 
is "?imax ' = (™85km)© ■ (-Hisskm)©- We determine 

/i max by scaling this equation. 

Values of h max for the different models are given in 
Table [3 Similarly to neutral elements absorbing light be- 
low /i max , it is likely that ionized elements are absorbing 
light above this limit, though we do not include this effect 
in the model. 

2.2.4. Presence of clouds 

In the wavelength range of interest, the surface of Venus 
is almost completely hidden by clouds. Therefore, it seems 
reasonable to model these types of clouds to a first order 
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Fig. 3. Atmospheric profiles, Bl. The legend is identical to 
that in Fig. [5] The temperature profile and mixing ratios 
are that of Venus. The temperature is considered to be 
constant above 50 km. Carbon dioxide is barely visible 
because it is by far the major constituent so its line is 
superimposed with that of the total density. 
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Fig. 4. Atmospheric profiles, CI. Same legend as in Fig. [3 
and Fig. The temperature profile follows a dry adiabat 
in the first 10 km of the atmosphere, until the point where 
e > e sat . Next, it follows a steeper saturated adiabat up 
to 20 km high. The temperature gradient is arbitrarily set 
to be isothermal above this point. The cloud top (thick- 
est line) is one scale height above the higher point where 
e > e sat . For reasons detailed in the text (see Sect. 12. 2. 5)1 , 
this point corresponds to the level where the temperature 
gradient becomes isothermal. 
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approximation by assuming that they act as an optically 
thick layer at a given altitude. As a result, clouds effec- 
tively increase the apparent radius of the planet and the 
transiting spectrum gives information only about atmo- 
spheric components existing above the cloud layer. The 
top of the cloud layer is a free parameter for N2/O2- and 
C02-rich atmospheres (set to 10 and 30 km, taken from 
the Earth and Venus, respectively). We treat the case of 
the N2/H20-rich atmosphere separately because H2O is a 
highly condensable species. 

2.2.5. Composition, vertical structure and location of 
the clouds in a IS^/H^O-rich atmosphere 

The temperature gradient of an atmosphere containing 
non-negligible amount of condensable species, like H2O, 
significantly departs from the case where no condensation 
occurs. A correct estimation of the temperature profile is 
crucial to determine the scale height, hence the detectabil- 
ity of that atmosphere. In an H^O-rich atmosphere, the 
evolution of the adiabatic temperature gradient is driven 
by the ratio of the partial pressure of water vapor, e, to 
the saturating vapor pressure, e sat . This ratio should also 
determine the levels at which the water vapor is in excess 
in the air and condenses (for e/e sa t > 1), i.e. the levels 
where clouds may form. 

Our initial conditions at the z = level (z°) are the 
temperature T° and pressure p . With these quantities we 
can estimate e sat , which depends only on the temperature, 
using the Clausius-Clapeyron equation: 

Hu 2 oLv f 1 1 

It* 



e sa ,t(T) =p*exp 



Af A k 



T 



(5) 



where p* and T* are the reference pressure (1.013 • 10 5 Pa) 
and temperature (373 K), /iH 2 o is the molar mass of wa- 
ter and L v is the latent heat of vaporization for water 
(2.26 • 10 10 ergg -1 ). Assuming that the planet is covered 
with liquid water (e.g., an ocean-planet; see Leger et al. 
2004) and that T° is 'tropical' (e.g. 340 K), the humid- 
ity at the surface is high so that the value of e° must be 
an important fraction of e sat (T°). We set e° to half the 
value of e sat (T°). The volume mixing ratio of water can 
be expressed as 1h 2 o = e/p, and we can calculate it at the 
surface of the planet. The atmosphere of an ocean-planet 
may also contain a significant quantity of CO2. We arbi- 
trarily set this quantity constant to Yco 2 =0.1 (Leger et 
al. 2004; Ehrenreich et al. 2005b). Molecular nitrogen is 
the major constituent of the atmosphere of the Earth and 
the second more abundant species in the atmosphere of 
Venus, and therefore we chose to include it to complete 
the chemical composition of this atmosphere. The mix- 
ing ratio of N 2 was set to be Y N2 = 1 — Y C q 2 — 1h 2 o 
at any level. Assuming the atmosphere contains only N 2 , 
H2O and CO2, we can obtain the mean molar mass of the 
atmospheric gas (/i° = ^ 3^° Mi) and that of the dry at- 
mospheric gas (/x^ = /i° — ^H 2 o/ i H 2 o) ; the mean specific 
heat of dry air (C° = ^ C Pi Yf/j,i/fj^) and the scale height 
H (all at the level z°). 



For the z 3+1 level, we need to evaluate the temperature 
gradient between z° and z J+ . There are two cases (Triplet 
& Roche 1986): 

— e J < e J sat ; in this case the temperature follows a dry 
adiabatic gradient, 



AT drj 



-9 
Cp 

in this case the gradient is saturated, 



(G) 



AT sat — ATdr\ 



(7) 



where r J sat 



(MH 2 oe4t)/[/4(y ~ e sat)] is the mixing- 
ratio of saturated air, R 3 dry = N A k / Mdry an( i ^h 2 o = 
N A k / ^h 2 o are the specific constant of dry air at the 
level z 3 and water (respectively). 

If z 3+1 < 20 km, we select the appropriate gradient ac- 
cordingly to the value of e/e sat , and get the value of the 
temperature T J+1 . Above 20 km, we assume the temper- 
ature profile becomes isothermal (T 3+1 = T 3 ). 

The assumption of an isothermal atmosphere, already 
discussed in Sect. 12.2.21 is somewhat arbitrary but is mo- 
tivated by an analogy with the atmosphere of the Earth, 
where the temperature gradient becomes positive from 
about 20 to 50 km. Taking an isothermal temperature gra- 
dient will conservatively mimic the presence of a strato- 
sphere. However, it has important consequences since it 
allows H 2 to be significantly present above the cloud 
top. In fact, above 20 km, the temperature stops decreas- 
ing, preventing condensation from occurring (the satura- 
tion vapor pressure depends only on temperature). Our 
assumption consequently fixes the height of the cloud deck 
to the point where the temperature profile is isothermal 
(actually, one scale height above that point). If we set 
this point higher, we would increase the amount of clouds 
hence reducing the detectable portion of atmosphere. In 
addition, the cloud formation would certainly take the cor- 
responding latent heat of condensation out of the atmo- 
spheric gas, contributing, as a consequence, to cool the 
atmosphere at the level of the cloud layer. 

We calculate H j+1 , p 3+1 = p> ■ exp(-z 3+1 /H 3+1 ), 
e 3 S!lX (from Eq. 03) and either e 3+1 = e 3 • 
exp [(z 3 ' — z 3+1 ) / H 3+1 )\ , if the atmosphere is not 
saturated or e 3+1 — eit. , if the atmosphere is saturated. 
We finally find all Y? , /u^L and and then iterate 

the process for all atmospheric levels. 

The higher and the lower pressure levels where e = e sat 
indicate respectively the bottom and the top of the region 
where clouds are forming. We assume the cloud layer does 
not extend over one scale height above the top of the cloud 
forming region. However, we can still have e < e sa t higher 
in the atmosphere, and thus H2O can be present above 
the clouds. 
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2.3. Description of atmospheric absorptions 

2.3.1. Chemical species 

We used the program LBLABC (Meadows & Crisp 1996), a 
line-by-line model that generates monochromatic gas ab- 
sorption coefficients from molecular line lists, for each of 
the gases, except ozone, present in the atmosphere. The 
line lists are extracted from the HITRAN 2000 databank 
(Rothman et al. 2003). We calculated the absorption co- 
efficients for O2, H 2 and CO2 in our wavelength range 
we (i.e., from 200 to 2 000 nm). 

The absorption coefficients relative to these species de- 
pend on pressure and temperature. We verified that those 
variations do not impact significantly on the results ob- 
tained (see Sect. 0} and we decided to use the absorp- 
tion coefficients calculated at the pressure and temper- 
ature of the cloud layer, i.e., 10 km in models Al, A2 
& A3, 30 km in models Bl, B2 & B3 and from 25 to 
70 km in models CI to C3. We then assumed these ab- 
sorption coefficients to be constant along the z-axis. This 
is a fairly good approximation since molecules at that at- 
mospheric level contribute more substantially to the trans- 
mitted spectrum than molecules at the bottom of the at- 
mosphere. Absorption coefficients for H2O, CO2, O3 and 
O2 are compared in Fig. 

The spectrum of O3 is unavailable in HITRAN at 
wavelengths lower than 2.4 /im. However it has strong 
absorption in the Hartley (200-350 nm) and Chappuis 
(400-750 nm) bands. Thus we took the photo- absorption 
cross-sections, a (in cm 2 ), from the GEIS A/cross-sectional 
databank (Jacquinet-Husson et al. 1999) and converted 
them into absorption coefficients, A (in cm 2 g _1 ), such as 
A = (tA/a/a<, where \i is the molar mass of the component. 

As shown in Fig. El the pressure and the temperature 
variations do not have a significant influence over the cross 
sections/absorption coefficients of O3. We therefore used 
the values given for p = 1 atm 3 and T = 300 K, and set 
them constant along the z-axis. 

2.3.2. Rayleigh diffusion 

Light is scattered toward short wavelengths by atmo- 
spheric molecules whose dimensions are comparable to 
A. Rayleigh diffusion could be an important indicator of 
the most abundant atmospheric species. Molecular nitro- 
gen, for instance, does not present any noticeable spectro- 
scopic lines between 0.2 and 2 /jm. With a transit obser- 
vation, the presence of a gas without spectroscopic lines 
like nitrogen in the Earth atmosphere can be indirectly 
inferred from the wavelength-dependance of the spectrum 
ratio continuum. Since Rayleigh scattering cross section of 
CO2 is high, Venus-like atmospheric signatures should also 
present an important Rayleigh scattering contribution. 

We have therefore estimated these different contribu- 
tions. The Rayleigh scattering cross section, <jr, can be 




Fig. 5. Absorption coefficients of atmospheric absorbers 
(in cm 2 g _1 ), as a function of the wavelength. The photo- 
absorption coefficients corresponding to H 2 0, O2, O3 
and CO2 (solid lines) are plotted against their respective 
Rayleigh scattering coefficient (dotted line), except O3, 
plotted against the Rayleigh scattering coefficient of N2 . 




270 

wavelength (nm) 

Fig. 6. Dependence of the absorption coefficient of O3 on 
pressure and temperature. For clarity, each line has been 
shifted down by 5 TO 4 cm 2 g _1 with respect to the previous 
one. 

expressed in cgs units as: (Bates 1984; Naus & Ubachs 
1999; Sneep & Ubachs 2004) 



vr{v) = 



r{uf 



(8) 



1 atm = 1 013 hPa. 



where v = 1/A, n is the number density (cm -3 ) and r is 
the refractive index of the gas. The total Rayleigh scat- 
tering includes weighted contributions from N2, O2, CO2 
and H2O (i.e., or = Y^i^i a Ri)i an d so we need all the 
corresponding refractive indexes. These are found in Bates 
(1984) and Sneep & Ubachs (2004) for N 2 , 2 and C0 2 . 4 

4 We noted a typographical error in the CO2 refractive in- 
dex formula (Eq. 13) in Sneep & Ubachs (2004): in order to 
yield the correct values, results from this expression should be 
divided by 10 3 (M. Sneep, personal communication). 
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The refractive index for H2O comes from Schiebener et al. 
(1990). Tests have proved the different refractive indexes 
do not significantly change with temperature and pressure. 
We have therefore calculated the indexes for standard con- 
ditions ( 15° C and 1013 hPa). 

2.3.3. Refraction 

Depending on the wavelength, the refraction may bring 
into the line of sight rays coming from different parts of 
the star. To quantify the importance of that effect, we 
calculate the maximum deviation, A9, due to the wave- 
length dependence of the refraction index, using the for- 
mula given by Seager & Sasselov (2000) and the refractive 
index at the surface (h = 0) between 0.2 and 2 /jm. We ob- 
tain A9 w 0.3'. This represents about 1.5%, 1% and 0.5% 
of the angular diameter of the star (F-, G- and K-type 
star, respectively) as seen from the planet. We can there- 
fore consider this effect negligible as long as there are no 
important variations of the stellar flux on scales lower than 
the surface corresponding to these numbers. 

2.4. Choice of test models 

We chose 9 cases, divided into 3 categories: 1 R ffi -planets 
(models Al, Bl and CI), 0.5 R e -planets (A2, B2 and C2) 
and 2 R^-planets (A3, B3 and C3). The parameters for 
each model are summarized in Table [5] For theses ranges 
of planetary radii, the depth of the occultation by the 
tested planets will differ by a factor of ~16 at most during 
their transit. Notice that a better detection of the tran- 
sit itself does not always imply a better detection for the 
atmosphere of the transiting planet. On the contrary, in 
some cases, the fainter the transit is, the more detectable 
the atmosphere will be! In any case, we naturally need to 
secure the detection of the planet itself before looking for 
an atmosphere. 

The choice of studying planets with a variety of sizes 
gives us the possibility to explore a large range of planet 
characteristics, in mass, radius and density. The Earth 
density is 5.5 gcm~ 3 . A planet having the internal com- 
position of the Earth and twice its radius would weigh 
~10 times more, while a planet half large would weigh 
~10 times less (Sotin et al. 2005). That gives densities 
of 6.1 and 4 gem -3 , respectively. We thus have 3 cases, 
each of which can be coupled with a plausible atmosphere. 
We chose a N2/02-rich atmosphere (similar to that of 
the Earth) for models Al, A2 and A3, and a Cytherean 
(i.e., Venus-like 5 ) COvrich atmosphere for models Bl, B2 
and B3. 

Note that the atmospheric pressure profiles are scaled 
from the 1 R ffi cases (Al and Bl) to the 0.5 and 2 R ffi 
models. In doing so, we did not include any species that 
showed a peak of concentration in altitude, such as the 

5 Cythera (Kv8i]pa) is an Ionian island where, according to 
the Greek mythology, the goddess Aphrodite/ Venus first set 
foot. See http://en.wikipedia.org/wiki/Cytherean. 
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O3 layer in model Al. In fact, the O3 peak does not de- 
pend only on the hydrostatic equilibrium, but also on the 
photochemical equilibrium at the tropopause of the Earth. 
For that reason O3 is absent in models A2 and A3. 

Leger et al. (2004) suggested the existence of 'ocean- 
planets', whose internal content in volatiles (H2O) might 
be as high as 50% in mass. Such planets would be much 
less dense than telluric ones. We are particularly inter- 
ested in those ocean-planets since the lower the density of 
the planet is, the higher the atmosphere extends above the 
surface. These objects could have densities of 1.8, 2.8 and 
4.1 gem -3 for radii of 0.5, 1 and 2 R ffi (Sotin et al. 2005), 
which are relatively small, but reasonable if compared 
with Titan's density (1.88 gem -3 ). The huge quantity of 
water on the surface of an ocean-planet could produce a 
substantial amount of water vapor in their atmosphere, 
if the temperature is high enough. A non-negligible con- 
centration of CO2 might be present as well in those at- 
mospheres (Ehrenreich et al. 2005b). Using this informa- 
tion on ocean-planets, we can simulate three extra cases, 
namely CI, C2 and C3 (Table 0). 

2.5. Choice of different stellar types 

In this work, we consider planets orbiting in the habitable 
zone (HZ) of their parent star. Our atmospheric models 
are not in fact a good description for planets orbiting too 
close to their parent star. For instance, the heating of the 
atmosphere by an extremely close star could trigger effects 
like evaporation, invalidating the hydrostatic equilibrium 
we assumed (see, for instance, Lecavelier des Etangs et al. 
2004; Tian et al. 2005). The reduced semi-major axis a r 
of the orbit of all planets we have considered is defined as: 

a r = a-(LjL Q )-°- 5 . (9) 

We set a r = 1 astronomical unit (AU) , so that the planet 
is in the HZ of its star. 

Here we focus on Earth-size planets orbiting around 
different main sequence stars, such as K-, G- and F- 
type stars, since the repartition of stellar photons in the 
spectrum is different from one spectral type to another. 
Planets in the HZ of K, G and F stars, with a r — 1 AU, 
should have a real semi-major axis of 0.5, 1 and 2 AU, 
respectively. 

3. Signal-to-noise ratio for ideal observations 

Prior to the atmospheres, we need to detect the planets 
themselves with a dedicated survey, as the one proposed 
by Catala et al. (2005). The transmission spectroscopy we 
theoretically study here require the use of a large space 
telescope. Hence, we need to quantify the S/N of such 
observations to determine the detectability of the atmo- 
spheric signatures for a transiting Earth-size exoplanet. 
The S /N will depend on both instrumental and astrophys- 
ical parameters. 
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3.1. Instrumental requirements 

The first relevant parameter relative to the instrumenta- 
tion is the effective area of the telescope collecting mirror, 
S, which can be expressed as S — (eD) 2 ir/A. The coeffi- 
cient e 2 accounts for the instrumental efficiency and eD is 
thus the 'effective diameter' of the mirror. Up to present, 
all exoplanetary atmospheric signatures have been de- 
tected by the Space Telescope Imaging Spectrograph 
(STIS) on board the Hubble Space Telescope (HST). This 
instrument, now no longer operative, was very versatile 6 
and consequently not planned to have high efficiency. It 
had a throughput e 2 w 2% from 200 to 300 nm, and 
e 2 f» 10% from 350 to 1 000 nm. As the majority of pho- 
tons we are interested in is available in the range from 350 
to 1 000 nm, we reasonably assume that a modern spectro- 
graph has a mean e 2 significantly greater than 10% from 
200 to 2 000 nm. Present day most efficient spectrographs 
have e 2 ss 25% in the visible, so it seems reasonable to 
imagine that next generation spectrographs, specifically 
designed to achieve high sensitivity observations, could 
have throughput of e 2 « 25%, or e = 50%. 

Another parameter linked to the instrument is the 
spectral resolution, TZ. In the following, 1Z will be assumed 
to be about 200, i.e. 10 nm-wide spectral feature can be 
resolved. 

Finally, it is legitimate to question the ability of the 
instrument detectors to discriminate the tenuous (~ 10~ 6 ) 
absorption features in the transmitted spectra of Earth- 
size planets. In a recent past, sodium was detected at a 
precision of 50 parts-per-million (ppm) on a line as thin 
as about 1 nm by Charbonneau et al. (2002) using STIS. 
According to our results (see Sect.QJ, some absorption fea- 
tures from Earth-size planet atmospheres show a ~1 ppm 
dimming over ~100 nm: the technological improvement re- 
quired to fill the gap should not be unachievable. Besides, 
since we deal with relative measurements - the in-transit 
signal being compared to the out-of-transit one - there is 
no need to have detectors with a perfect absolute calibra- 
tion. Only a highly stable response over periods of several 
hours is required. Nevertheless, instrumental precision re- 
mains a challenging issue whose proper assessment will 
require further, detailed studies. 



3.2. Physical constraints on the observation 

The number of photons detected as a function of wave- 
length depends on the spectral type of the star, while the 
total number of photons received in an exposure of dura- 
tion t depends on the apparent magnitude of the star, V. 
The stellar spectra FY =0 (X) are from p Capricorni (F2 iv), 
HD 154760 (G2 v) and HD 199 580 (K2iv) and are taken 
from the Bruzual-Persson-Gunn-Stryker (BPGS) spec- 
trophotometry atlas 7 . The fluxes (ergcm~ 2 s _1 A -1 ) are 
given at a null apparent magnitude, so we re-scaled them 




1000 
wavelength (nm) 



Fig. 7. Spectrum of a K2 (dashed line), G2 (solid line), 
and F2-type stars (dotted line) between 0.2 and 2 pm. 
The fluxes are scaled to an apparent magnitude V = 8. 

for any apparent magnitude V, F+ = FY~° ■ 10 -0-4V '. The 
three corresponding spectra are plotted for a default mag- 
nitude V = 8 in Fig. □ 

The stellar type determines the radius and the mass 
of the star, so the transit duration (and thus the maxi- 
mum time of exposure during the transit) is different de- 
pending on the star we consider. The transit duration is 
also a function of the semi-major axis of the planet orbit. 
Since we chose a constant reduced distance (o r = 1 AU) 
for all planetary models (see Sect. 12.51) . the duration of 
transit depends on the stellar luminosity as well. From 
Zombeck (1990), we obtain the radii of F and K stars rel- 
atively to that of the Sun, respectively Ry/Rq « 1.25 
and Rk/Rq ~ 0.75, the mass ratios, M F /M Q « 1.75 and 
Mk/Mq w 0.5, and the luminosity ratios, respectively 
Ly/Lq w 4 and Lk/Lq w 0.25. Using Eq.[§l the duration 
of the transit is: 



13tt R* 



M* 

M 



0.25 



(10) 



6 STIS was used for imagery, spectro-imagery, coronography 
and low and high resolution spectroscopy. 

7 Available on ftp.stsci.edu/cdbs/cdbs2/grid/bpgs/. 



where 137r/4 h is the mean transit duration of a planet 
at 1 AU across a G star averaged over all possible impact 
parameter of the transit. From Eq. 1101 we obtain mean 
transit durations of 7.6, 10.2 and 13.6 h for K-, G- and 
F-type star, respectively. In the following, we set t = r. 

Ideally, our observations are limited only by the stel- 
lar photon noise - the detection of sodium at a preci- 
sion of ~50 ppm in the atmosphere of HD 209 458b by 
Charbonneau et al. (2002) was in fact limited by the stel- 
lar photon noise. However, at the low signal levels we are 
searching for, the intrinsic stellar noise might need to be 
considered as well. Stellar activity, as well as convective 
motions will cause variations in both intensity and color 
in the target stars, on a large variety of timescales. The 
impact of stellar micro-variability on the detectability of 
photometric transits has been addressed by a number of 
studies (see, e.g., Moutou et al., 2005; Aigrain et al., 2004 
- especially their Fig. 8; Lanza et al., 2004), all point- 
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ing towards photometric variability levels in the range of 
~ 100-1 000 ppm for durations of a few days. This is to be 
compared to the strength and duration of the atmospheric 
signatures we want to look at: they are ~1 ppm variations 
lasting a few hours. While indeed the different time fre- 
quency and spectral content of these signatures versus the 
stellar noise will hopefully allow to discriminate the two, 
the impact of stellar micro-variability on such faint sig- 
nals is likely to be significant, and may limit the ability to 
detect an atmosphere in a transiting planet. For instance, 
Aigrain et al. (2004) suggested K stars are more adapted 
than G or F stars regarding to the detection of terrestrial 
planets versus stellar micro- variability. However, note that 
the observation of several transits for each planet consid- 
ered will confirm the signal detected in the first transit. 
For instance, at a r = 1 AU around a K star, a planet has 
a period of « 0.3 yr, allowing to schedule several tran- 
sit observations within a short period of time. Finally, the 
usual technique to detect a spectral signature from a tran- 
sit is to compare in-transit and out-of-transit observations 
(Vidal-Madjar et al. 2003, 2004). For all these reasons, we 
will assume in the following to be able to discriminate a 
transit signal from the stellar activity and consequently 
the photon-noise to be the limiting factor. Nevertheless, 
further and detailed analysis is certainly needed to quan- 
tify the effect of stellar micro- variability, as a function of 
the stellar type, but this is outside the scope of this paper. 

3.3. Calculation of the signal-to-noise ratio 

Now let tp+ be the maximum number of photons per 
element of resolution that can be received during r: 
y>* = F is (X) ■ A/ (hpc) ■ TZ ■ S ■ t, where hp is Planck's con- 
stant and c the speed of light. Some photons are blocked 
or absorbed by the planet, therefore the actual number of 
photons received during the transit is ip = + 5ft') per 
element of resolution. 

From the observations, it is possible to obtain Rp, an 
estimate of the radius of the transiting planet Rp (e.g., 
by using the integrated light curve or a fit to the observed 
spectrum ratio). This value corresponds to the flat spec- 
trum ratio (i.e., a planet without atmosphere) that best 
fits the data. The corresponding number of photons re- 
ceived during an observation per element of resolution is 

therefore expressed as: tp — ip+ 1 — (Rp/R+) 2 . 

The weighted difference between ip and tp can reveal 
the presence or the absence of a planetary atmosphere. 
We express the \ 2 of this difference over all the elements 
of resolution k as J^k [(V* — Vk) / a vtX Here, the uncer- 
tainty of the number of photons received is considered to 
be dominated by the stellar photon noise (see Sect. I3.2|l . 
that is a v = yftp. We thus have: 



X 




(hp/R* 



1 2 N 



(11) 



Given the \ 2 , the S/N can be directly calculated taking 
its square root. The best estimation can be obtained by 



minimizing the x 2 with respect to the radius Rp, i.e., 
dx 2 /dRp = 0. From this formula we can calculate the es- 
timated radius: 



R P = R+* 



E*[vW(i + «L)] ' 



(12) 



Once wc determine if an atmosphere is observable or 
not (depending on the S/N ratio), we can use a similar 
approach to quantify the detectability of the single atmo- 
spheric absorber contributing to the total signal ip. Let 
(pi = ip*(l + 3?'i) be the signal obtained by filtering the 
stellar light out of all atmospheric absorbers except the 
i th , and let {(pi) be its estimation. Here, Sft- is the spec- 
trum ratio calculated when the species i is not present in 
the atmosphere. Further, since (ipi) « ontpi, we can deduce 
the presence of absorber i in the atmosphere, by simply 
comparing the fit we made assuming its absence (oti<pi) 
with the measured signal (tp): 



X] 



E 



i + K 



(l + 5ft' fe )-a l (l + 5R' !fc )] 2 V (13) 



where 





^ fe (i + 34) 






v+k ( 


i + ^ fc ) 2 /(i 





(14) 



4. Results and discussion 

The results of our computations are displayed in Tables [3] 
fe^and plotted as spectrum ratios in Figs. 151 151 fclTUl 



4.1. Spectral features of interest 

Here we summarize the contributions of each atmospheric 
absorber to the spectrum ratio for various models. The 
spectral resolution of the plots presented here is 10 nm. 
The most prominent spectral signatures, when present, 
are those of O3 and H 2 0. Carbon dioxide is hard to dis- 
tinguish from H 2 bands and/or its own Raylcigh scatter- 
ing. Molecular oxygen transitions are too narrow to sig- 
nificantly contribute to the spectrum ratio. 

4.1.1. Ozone 

In the spectral domain studied here, the Hartley (200- 
350 nm) and Chappuis (420-830 nm) bands of O3 ap- 
pear to be the best indicators of an Earth-like atmosphere. 
These bands are large (respectively 150 and 600 nm) and 
lay at the blue edge of the spectrum, where spectral fea- 
tures from other species are missing. There is noticeably 
no contamination by H2O, and O2 strong transitions are 
narrow and could be easily separated. Ozone bands signif- 
icantly emerge from Rayleigh scattering and they corre- 
spond to very strong transitions, despite the small amount 
of O3 present in the model Al atmosphere (Y"o 3 < 10~ 5 ). 
When present, ozone is more detectable in an atmosphere 
similar to model A2. 
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4.1.2. Water 

The signature of H2O is visible in a transit spectrum only 
if H2O is substantially abundant above the clouds. This 
is not the case for models of Earth-like atmosphere like 
Al, A2 and A3. On the contrary, the models of the ocean- 
planets (CI, C2 and C3) show a major contribution from 
this molecule, in the form of four large bands that dom- 
inate the red part of the spectrum (at A ^ 950 nm). 
For these three cases, H 2 can be significantly abundant 
above the clouds. 

4.1.3. Carbon dioxide 

The lines of CO2 are about as strongly emerging from the 
'continuum' than the H2O ones, but are often overlapping 
with these lines. The transitions around 1 600 nm and the 
ones around 1 950 nm are the easiest to identify, other 
bands are not observable if water is present. Rayleigh 
scattering and photo-absorption cross sections of C0 2 are 
comparable at most wavelengths below 1.8 /mi (see Fig.QjJ), 
except for a few ~10-nm wide bands. In fact, the more 
CO2 is present in the atmosphere, the more opaque the 
atmosphere becomes. This implies it would be impossible 
for an observer on the surface of Venus to see the Sun. 
Carbon dioxide may be more detectable farther in the in- 
frared, hence making desirable further investigations up 
to 2.5 /im. 

4.1.4. Molecular oxygen 

Molecular oxygen does not appear in the plots: its bands 
at 620, 700, 760 and 1 260 nm are too thin to appear 
with only 10 nm resolution. Besides, its Rayleigh scatter- 
ing cross section almost completely masks its absorption 
features (see Fig. [S| so that n0 large bands of O2 can be 
used as an indicator of its presence. However, note that 
the presence of O3 indirectly indicates the presence of O2, 
as pointed out by Leger et al. (1993) and others. 

4.1.5. Rayleigh scattering 

When Hartley and Chappuis bands of O3 are absent (all 
cases but Al), the Rayleigh scattering signature is clearly 
visible in the blue part of the spectrum ratio. On one side 
it masks the presence of some transitions, like those of O2 
and some of CO2, but on the other side it can provide 
two important informations: (i) even if the spectral fea- 
tures cannot be distinguished because they are too thin 
or faint, the characteristic rising 'continuum' as A~ 4 for 
short wavelengths is a clear indication that the planet has 
an atmosphere, and (ii) it indirectly indicates the presence 
of the most abundant species of the atmosphere, such as 
CO2 and N2, even if N2 shows no spectral signature in 
the observed domain. As a consequence, Rayleigh scatter- 
ing can be considered a way to detect N2, provided clouds 
and/or aerosols do not in turn mask the Rayleigh scatter- 
ing signature. 



To summarize, it is possible to detect the presence of 
the atmosphere of a transiting cxoplanct thanks to the 
Rayleigh scattering, whatever the composition of the at- 
mosphere is. Moreover, it is theoretically possible to dis- 
criminate between an 02-rich atmosphere, where O3 is ex- 
pected to be present (Leger et al. 1993; Sagan et al. 1993) 
and a H2 0-rich atmosphere, as the O3 lifetime is supposed 
to be extremely brief in a water-rich environment. In other 
words, we should be able to distinguish telluric Earth-like 
planets with low volatile content from volatile-rich plan- 
ets. On the other hand, high spectral resolution is needed 
to discriminate between H2 0-rich planets and Cytherean 
worlds (Bl, B2, B3). 

4.2. Parameters influencing the signal-to-noise ratio 
4.2.1. Influence of the star 

From Table 13 it is clear that the best targets are K-type 
stars, rather than G- or F-type stars, the former allowing 
much better S/N than the latter. Two factors are deter- 
mining the role of the star in the capabilities of detect- 
ing an exoplanet atmosphere: (i) The size i?* of the star, 
which directly influences the S/N (see Eq. Illf) and the 
duration of transit (Eq. HOfl and (ii) the semi-major axis 
of the planet's orbit, which influences both the duration 
of transit and the probability to observe the transit from 
Earth (see below). These factors can explain the discrep- 
ancies between the S/N values obtained for different kind 
of stars in Table El 

The probability, a, that a planet transiting its par- 
ent star might be seen from the Earth is defined as 
a = _P{transit} = R±/a, with i?* being the radius of 
the star and a the semi-major axis of the planet's orbit. 
This probability is about 10% for 'hot Jupiters', while it 
is 0.3%, 0.5% and 0.7% for planets orbiting in the HZ of 
a F, G or K star, respectively. 

In addition, K stars are more numerous than other 
types of stars. From the CDS database, we find there is 
approximatively a total of 10 OOOTO ' 6 ^^ 8 -' main sequence 
stars brighter than a given magnitude V on the whole 
sky. 8 About 3/5 of these are K type stars, against only 
1/10 for G stars. Let us now define f3 to be the number 
of planet(s) per star, and 7 to be the fraction of the sky 
that is considered for a transit detection survey (in other 
words the efficiency of surveys to find the targets). We list 
in Table0]the number of potential targets for each model. 
This number, N, corresponds to the number of targets 
detected with a 10-m telescope mirror effective size and 
with a S/N greater than or equal to 5. It is given by: 

AT AT Q /S/Ny_g m \ 

Ns/N>5, £ fl=iom = No-a-fi-^- 1 - I , (15) 

where iV is about 6 000, 1000 and 3 000 for K, G and F 
stars respectively, i.e. the number of stars with magnitude 

8 We consider mostly bright stars, for which the distribution 
is essentially isotropic. 
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> 8, and S/Nv=8,eD=i0m is the expected S/N ratio com- 
puted for a given atmosphere of a planet orbiting a V = 8 
star with a telescope having a mirror effective area of 10 m 
(this value is given in the last column of Tabled . Since no 
Earth-size planet has been discovered so far, we have no 
real estimate of (3. In the following, when it is not a free 
parameter we consider (3 = l 9 . Catala et al. (2005) pro- 
pose a 30° x 30° survey dedicated to find planets around 
< ll th -magnitude stars, i.e., 7 « 2-3% for such a project. 

Let be iV s / N eD , the number of potential targets reach- 
ing a minimum S/N ratio for a given mirror effective size 
eD, which scales from the value calculated using Eq. 
^ r s/N>5,e/3=i0m, in the following way: 

/S/N\ -3 / eD \ 3 

Ns/N, eD =iVs/N>5, e£>=10m " ( — I ' I 10~m / ^ 

The values obtained for atmospheric detection are 
strongly in favor of a small, late type star. Note that this 
is also true for the detection of the planetary transit as 
well. 

4.2.2. Effect of the atmospheric temperature gradient 

The thick CO2 Venus-like atmospheres (Bl, B2 and B3, 
see Table |3 & EJ are more difficult to detect than other 
cases. Even if we set the top of the clouds at 10 km 
height, the detection remains more challenging than for 
model Al. That is somewhat surprising, partly because 
CO2 has strong transitions, particularly in the near in- 
frared, and partly because of the larger scale height at the 
surface of the planet (14.3 km for model Bl, 8.8 km for 
model Al). As a consequence, the atmosphere in model Bl 
should have a larger vertical extent than in model Al. 
In reality, the difficulty to characterize the atmospheres 
of models Bl, B2 and B3 is related to the temperature 
profiles we chose (see Fig. |3 and OJ: At 50 km of al- 
titude, the temperature of model Bl is roughly 60 K 
colder than that of model Al. This model in fact ben- 
efits from the positive stratospheric temperature gradi- 
ent of the Earth. Moreover, the atmosphere for model Bl 
(Mbi = 43 gmol -1 ) is heavier than the one for model Al 
(a*ai = 29 gmol -1 ). Therefore, at high altitude, the scale 
height is larger in model Al than in model Bl (respec- 
tively 7.6 km and 3.9 km at an altitude of 50 km). 

4.2.3. Effect of atmospheric pressure 

Note that the thickness of the atmosphere in model Bl 
is almost half the one in Al, despite the intense surface 
pressure (100 atm), which should help to increase the up- 
per level of the atmosphere, limited by the UV photo- 
dissociation (h max ). The exponential decrease of pressure 
prevents, in fact, po to play a key role: in order to coun- 
terbalance the effect of the negative temperature gradient, 

9 Actually, j3 — 2 in the Solar System because there are two 
Earth-size planets with atmospheres, namely Venus and the 
Earth. 



the surface pressure should have been > 10 6 atm to obtain 
absorptions similar to the case of the Earth (model Al). 

4.2.4. Effect of the planet gravity and density 

The atmospheric absorption is, at a first order, propor- 
tional to H Rp. At a given temperature and for a given at- 
mospheric composition, the scale height H is proportional 
to the inverse of the gravity acceleration, g —1 , or equiv- 
alently to R P /Mp, where Mp is the mass of the planet. 
As a result, the absorption is expected to be roughly in- 
versely proportional to the bulk density of the planet, pp, 
independently of the planet size. 

This effect is illustrated by the following examples: 
models CI, C2 and C3 all benefit from very extended 
atmospheres, given the weak value of g in the three 
cases. For a planet as dense as the Earth (such that 
9ci — <?ai), the results for the N2/H20-rich atmosphere 
in models C are close to the ones obtained for models A. 
Both models C and A, present typical spectral features. 
In model Al, ozone, for which the concentration peaks at 
the tropopause, gives a prominent signature in the blue 
edge of the spectral domain (the Hartley and Chappuis 
bands, as seen in Fig. top panel). On the contrary, the 
saturated atmosphere of model CI, which sustain H2O up 
to high altitudes, yields strong bands around 0.14 and 
0.19 /jm (Fig. |S| bottom panel). The role played by g 
can be better understood by comparing model A3 or B3 
(g = 24.5 ms" 2 ) to model A2 or B2 (g = 3.9 ms -2 ), and 
model C3 (g = 14.7 ms^) to model C2 (g = 2 ms -2 ). 
Using absorption spectroscopy, it is clear that the atmo- 
spheres of small and light planets (i.e., with low surface 
gravity) are easier to detect than the ones of large and 
dense planets (i.e., with high surface gravity). 

Small and light exoplanets, however, may not be able 
to retain a thick atmosphere. In fact, high thermal agita- 
tion of atmospheric atoms causes particles to have a ve- 
locity in the tail of the Maxwellian distribution allowing 
them to escape into space (i.e., Jean's escape). It is there- 
fore questionable if planets of the size of Titan can have a 
dense atmosphere at 1 AU from their star. Models A2, B2 
and C2 enter that category. This problem concerns both 
small planets and giant exoplanet satellites. 

According to Williams et al. (1997), a planet having 
the density of Mars could retain N and O over more than 
4.5 Gyr if it has a mass greater than 0.07 M$. Model 
planets A2 and B2 have masses of 0.1 M ffi and a den- 
sity equivalent to that of Mars («4 gem -3 ) so they would 
be able to retain an atmosphere (though they may not 
be able to have a 1 atm atmosphere, as for Mars). The 
ocean-planet model C2 has a mass of 0.05 M ffi for a den- 
sity of 2.8 gem -3 , and according to Williams et al. (1997), 
its atmosphere should consequently escape. However, al- 
though at 1 AU from the star, such a planet also has a 
huge reservoir of volatile elements. This reservoir should 
help to 'refill' the escaping atmosphere. 
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Note that an hydrodynamically escaping atmosphere 
should be easier to detect than a stable one, since it can 
bring heavier elements into the hot upper atmosphere. 
This effect is illustrated by the absorptions seen by Vidal- 
Madjar et al. (2003, 2004) in the spectrum of HD 209458, 
which originate in its transiting giant planet hydrodynam- 
ically escaping atmosphere. A model of an 'escaping ocean' 
is studied by Jura (2004). This process would give interest- 
ing absorption signatures in the H2O bands from the lower 
atmosphere and in the signatures of the photo-dissociation 
products of H 2 from the upper atmosphere (such as an 
absorption of Lyman a photons by the hydrogen atom). 
See detailed discussion in Jura (2004). 



5. Conclusion 

The vertical extent of the atmosphere is of extreme im- 
portance as concerns the detectability of a remote atmo- 
sphere by absorption spectroscopy. This tends to favor less 
dense objects, like giant exoplanet satellites (as would be 
an 'exo-Titan') or volatile-rich planets (as ocean-planets, 
theoretically possible but not observed yet). Cytherean 
atmospheres are the most challenging to detect. Surface 
parameters, such as surface pressure and temperature, are 
not crucial. A temperature gradient that becomes positive 
at few tens of kilometers height (for instance owing to pho- 
tochemistry) might help the detection. Our results show 
that late-type stars are better for detecting and character- 
izing the atmospheres of planets in transit, since they are 
smaller, more numerous and present a better probability 
of being transited by a planet. 

The strongest signatures of the atmosphere of a tran- 
siting Earth-size planet could be those of H2O (6 ppm in 
the case of hypothetical ocean-planets), O3 (~l-2 ppm) 
and CO2 (1 ppm), considering our spectral study from 
the UV to the NIR (i.e., from 0.2 to 2 (im). The pres- 
ence of an atmosphere around hundreds of hypothetical 
'ocean-planets' (models C) could be detected with a 10- 
20 m telescope. The atmospheres of tens of giant exoplanet 
satellites (model A2) could be in the range of a 20-30 m in- 
strument. A 30-40 m telescope would be required to probe 
Earth-like atmospheres around Earth-like planets (model 
Al). These numbers suppose that Earth-size planets are 
frequent and are efficiently detected by surveys. 

Finally, planets with an extended upper atmosphere, 
like the ones described by Jura (2004), hosting an 'evapo- 
rating ocean', or the planets in an 'hydrodynamical blow- 
off state', are the natural link between the planets we have 
modelled here and the observed 'hot Jupiters'. 
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Fig. 8. Spectrum ratios for models Al (a), Bl (b) and CI 
(c) . The spectrum ratios have been respectively shifted by 
the values in parenthesis so that the absorption by the 
'solid disk' of the planet is ppm. In the case of models 
with clouds, the 'solid disk' is artificially increased by the 
cloud layer. The dashed line indicates the best fit estima- 
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Table 2. Summary of test models. 
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Table 3. Summary of results: signal-to-noise ratios obtainable with a telescope mirror effective size of eD = 10 m 
pointing at a V = 8 star. To get the S/N ratios for a different effective size eD, exposure time during transit, t, 
and/or apparent magnitude of the star, V, the result scales with (eD/10 m) • (t/r) - 5 ■ lfj -0 - 2 ^ -8 ) where r is defined 
by Eq. The S/N by species are calculated for the models with clouds. 



References 

Aigrain, S., Favata, F., & Gilmore, G. 2004, A&A, 414, 1139 
Bates, D. R. 1984, P&SS, 32, 785 
Brown, T. M. 2001, ApJ, 553, 1006 

Catala, C., Aerts, C., Aigrain, S., et al. 2005, Proc. 39 th ESLAB 
Symposium, Noordwijk, 19-21 April 2005, Favata F., & 
Gimenez, A., eds. 

Charbonneau, D., Brown, T. M., Latham, D. W., et al. 2000, 
ApJ, 529, 45 

Charbonnneau, D., Brown, T. M., Noyes, R. W., & Gilliland, 
R. L. 2002, ApJ, 568, 377 



Dubois, V. 2002, PhD. thesis, Universite de Nantes - 

Observatoire de Paris 
Elliot, J. L., & Olkin, C. B. 1996, ARE&PS, 24, 89 
Ehrenreich, D., Despois, D., Leger, A., et al. 2005, in prepara- 
tion 

Fishbein, E., Farmer, C. B., Granger, S. L., et al. 2003, IEEE 
Transactions on Geoscience and Remote Sensing, 41, 2 

Green, D., Matthews, J., Seager, S., & Kuschnig, R. 2003, ApJ, 
597, 590 

Hubbard, W. B., Fortney, J. J., Lunine, J. I., et al. 2001, ApJ, 
560, 413 



14 



D. Ehrenreich et al.: Spectrum of Earth-size transiting planets 



Model 


Description 


Star 


Mirror 


Limiting 


Number 




Number of tar 


gets 








off. size (m) 


magnitude 


of stars 


for models w/ clouds 








V c ^/b/N>5, V— 8 

w/ clouds 


( Vljim) S /N>5, — 10 m 

w/ clouds 




0-7=1 
D = 20 rr 


(A0s/n>5, e = 
0-7 = 3% 
D = 30 m 


50% 
P ■ 7 = 10% 
D = 30 m 








(a) 


(b) 


(c) 




(d) 








K 


14 


7.22 


2 042 


14 


1 


4 


Al 


(«)Earth 


G 


22 


6.31 


96 


< 1 (0.4) 


< 1 


< 1 (0.1) 






F 


29 


5.66 


118 


< 1 (0.3) 


< 1 


< 1 (0.1) 






K 


21 


6.31 


580 


4 


< 1 (0.4) 


1 


Bl 


(~)Venus 


G 


43 


4.90 


13 


< 1 


< 1 


< 1 






F 


68 


3.73 


8 


< 1 


< 1 


< 1 




medium 


K 


1.3 


12.5 


> 3 • 10 b 


19 602 


1984 


6 615 


CI 


ocean- 


G 


2.5 


11.0 


63 095 


321 


32 


108 




planet 


F 


3.9 


10.1 


54 591 


157 


15 


52 






K 


8 


8.50 


11 971 


84 


8 


28 


A2 


small Earth 


G 


13 


7.51 


508 


2 


< 1 (0.2) 


< 1 (0.6) 






F 


17 


6.90 


656 


1 


< 1 (0.1) 


< 1 (0.3) 






K 


15 


7.10 


1 730 


12 


1 


4 


B2 


small Venus 


G 


31 


5.52 


32 


< 1(0.1) 


< 1 


< 1 






F 


50 


4.50 


23 


< 1 


•C 1 


< 1 




small 


K 


1.1 


12.8 


> 4 ■ 10 b 


33 569 


3 398 


11329 


C2 


ocean- 


G 


2.0 


11.5 


125 892 


600 


60 


202 




planet 


F 


3.1 


10.5 


94 868 


307 


31 


103 






K 


45 


4.71 


63 


< 1 (0.4) 


< 1 


< 1 (0.1) 


A3 


super- Earth 


G 


86 


3.39 


1 


< 1 


< 1 


< 1 






F 


121 


2.51 


1 


< 1 


< 1 


< 1 






K 


> 10 3 







< 1 


< 1 


< 1 


B3 


super- Venus 


G 


> 10 3 







< 1 


< 1 


< 1 






F 


> 10 3 







< 1 


< 1 


< 1 




big 


K 


4.0 


10.1 


109 182 


682 


69 


230 


C3 


ocean- 


G 


7.7 


8.57 


2197 


10 


1 


3 




planet 


F 


13 


7.57 


1656 


4 


< 1 (0.4) 


1 



Table 4. Summary of results: mirror effective size and number of targets. 

(a) Effective size (e-D)s/N>5. v=8 °f the telescope mirror required to obtain S/N = 5 for a V = 8 star, based on the 
numbers displayed for the models with clouds (see Table 0). 

(b) The limiting magnitude at which the number of targets in the last column is given. This can be expressed as 



(VLim)s/N> 



5, eD=Wm 



8. 



5 • log 10 [(s/N y=8 eD=10 j /5 ■ (eD)/10 m 

(c) Total number of given spectral-type stars brighter than the limiting magnitude. 

(d) Number of potential targets calculated with Eq. using the S /N value of the models with clouds and assuming 
various f3 ■ 7 values. The coefficients (3 and 7 are defined in the text. When the number of potential targets is slightly 
less than 1, the value is given between parenthesis. Use Ea. llfil to scale the value displayed in the column to any mirror 
effective size eD and minimum S/N. 
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Fig. 9. Spectrum ratios for models A2 (a), B2 (b) and C2 
(c). The 'saturation effect' in H2O lines, for model C2, 
is a consequence of the atmosphere being optically thick 
at the upper atmospheric level, h max . In fact, if one con- 
sider there is no more water above this level due to photo- 
dissociation (see Sect. I2.2~3"l) . such transmitted spectrum 
plots allow to determine the level where H2O photo- 
dissociation occurs in an exoplanet atmosphere. 
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(c). 
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